Manifestation of superfluidity in an evolving Bose-condensed gas 
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We study the generation of excitations due to an "impurity" (static perturbation) placed into an 
oscillating Bose-condensed gas in the time-dependent trapping field, ft is shown that there are two 
regions for the position of the local perturbation. In the first region the condensate flows around 
the "impurity" without generation of excitations demonstrating superfluid properties. In the second 
region the creation of excitations occurs, at least within a limited time interval, revealing destruc- 
tion of superfluidity. The phenomenon can be studied by measuring the damping of condensate 
oscillations at different positions of the " impurity" . 
In studies of Bose-Einstein condensation in ultracold gases there has been outlined a new stage associated with 
investigating superfluidity (SF) in such systems. This problem is of a specific interest since it is a question of SF of 
a dilute gas. The first results in revealing a critical velocity have been obtained in At the same time in |^ the 
creation of a vortex is demonstrated in a sophisticated experiment with two condensates. Recently, direct creation of 
vortices has been observed in a rotating condensate 

An interesting possibility for studying SF may be realized in the observation of evolution of condensate accompa- 
nying time-dependent variations of the confining field. For the most interesting case when a parabolic potential keeps 
its shape, the scaling solutions found in |^-|^ allow us to describe the evolution of condensate for arbitrary variations 
of the trap frequencies. The analysis of solutions, in particular, given below, shows that during the evolution there 
are space-time regions where the velocity v{r,t) of condensate exceeds the local speed of sound c(r, t), i.e., the local 
critical velocity for the SF system. Simultaneously, there exists a spatial region where, on the contrary, v{f, t) < c(r, t) 
for any time. In the absence of external perturbations the presence of the regions with w(r, t) > c(r, t) does not result 
in itself in the irreversible processes and the oscillations of condensate do not decay at a temperature T = 0. 

However, provided an external perturbation is localized in the region where w(r, t) > c(r, t) even for a limited 
time interval, generation of collective excitations of the condensate takes place. Naturally, this results in the decay 
of the oscillations of the condensate. At the same time, if the perturbation is located in the region where always 
v{f,t) < c{f,t) , the gas in the presence of SF will flow adiabatically about the local perturbation not creating 
excitations and, therefore, without decay of the oscillations of the condensate. Thus, measuring the decay of oscillations 
of the condensate under changing localization of perturbation, we can display the SF properties of a Bose-condensed 
gas. 

As is found in [Q, for the isotropic 2D parabolic potential there is an exact scaling solution for arbitrary values 
of the gas density n and initial trap frequency loq and an arbitrary dependence i^(t). Moreover, the energy of the 
interparticle interaction and kinetic energy in this case change identically with an expanding gas, entailing similar 
variation of the main parameters. On the whole, this makes a choice of the trap geometry close to the cylindric 
one very attractive. In this case a static perturbation can be induced by the laser beam parallel to the longitudinal 
cylindric axis and limited in the transverse size. Shifting the beam axis in the radial direction, one can change the 
localization of perturbation. 

It is worthwhile to note that for this type of investigations, a trapped Bosc- condensed gas has an essential advantage 
in comparison with the *He quantum liquid. This is related with the isolation of a gas from any walls, the possibility 
to concentrate a majority of particles in the condensate, and with studying the appearance of irreversibility and the 
partial destruction of coherence at the conditions of the large-scale time-spatial evolution of condensate with the 
strong change of density 

Below we consider the evolution of the condensate in an isotropic 2D parabolic potential with the time-dependent 
frequency u (t) and determine the generation of excitations under external static perturbation. The generalization 
for the 3D case of the cylindric symmetry and extended perturbation unchanging along the longitudinal axis can be 
performed easily. 

Examining evolution of the 2D interacting Bose gas in the time-dependent confining potential, we can employ the 
general equation for the Heisenberg field operator of atoms 

i^i \ ' - -T^A + ^ * f, t) + C/o*+ r, i * f, i * f, t). 1 

at 2m 2 

The only simplification in the equation is an assumption of the local character of the interparticle interaction which 
allows us to present this interaction in the form U {ri — — UqS (ri — r2). 
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Let us introduce a scaling parameter b {t) and, correspondingly, spatial and time variables p = r^^Ot) , 
For the 2D case, let us represent the operator ^ (r, t) as 

«'(r,i) = ix(p,r)exp[z$(r,t)]. (2) 
Substituting (^) into (|l|) and using the results obtained in , we find 

"^^^'^-"-^mbdf 

Then the equation for operator x {p, t) reduces to 

dx 1 

= -^ApX+ -mtj^p^X + t^oX^XX- (4) 

The equation takes such form if we accept that b {t) and T{t) are determined by the equations 

t 

f^+co^t)b^^, rit) = jdt'/b\t'). (5) 



We put (— oo) = ujo in Eqs. ^ and (^. Equation (||) should be solved at the initial conditions of 6(— oo) = 1, 
6(-oo) = . 

It follows from Eq. that the problem in the p and r variables reduces to the solution of the equation for the 
static parabolic potential of the initial frequency wq- Finding solution and deriving b {t) and r [t) from (|^), we obtain 
a complete description for the space-time evolution at an arbitrary variation of uj (t) . It is essential that this refers, 
to the same extent, to both the condensate and the excited states of the system. 

Let us restrict our consideration with the case of zero temperature. Then the ground state in the static potential 
represents the condensate and in Eq. (|^) the operator x can be replaced by the macroscopic condensate wavefunction 
Xo (Pi ''') having a typical dependence on r 

Xo(p,T)=xo(p)e-^^^/^ (6) 
Here ^ is the starting chemical potential and xo {p) is the solution of the equation 

- ^^p^o + [-M + -j^mulp^ + C/oXolxo = 0. (7) 

To describe excitations in the system, we introduce operator x' (p, r) . Let us write the initial operator x 
X {p, t) = [xo (p) + x! {P: t)] e~^f^'^/^ and substitute it into (^. Then we find for the equation Hnearized in x' (p, t) 

= -^\X' + {^mu;^ p)x' + 2Uoxlx' + Uoxlx'^ ■ (8) 

This equation determines excitations in the coordinates of (p, t) and we can consider any problem in this reference 
frame. However, keeping in mind the solution of the problem on creating excitations under a static perturbation, it 
is more natural to investigate it in the laboratory reference frame. 

Let us rewrite Eq.(0) introducing new variables rand t. For this purpose, we use the relations 



dx' ,2<9x' , -^;dx' . ,2 , 



We assume that the density of a gas is sufficiently large and p S> huo, i.e., the Thomas-Fermi approximation is 
valid. Then we find the condensate density from Eq.(|^) (see (|^)) 

"o(r,t) = ^xg, = (9) 

where 
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R{t)^Rob{t), Ro = ^2^i/mujl. (10) 
As a result, in the laboratory reference frame Eq.(||) reads 

tn^^~—Ax'+gx' + 9x'+-invS/x'. (11) 

Here 

v = f^, (12) 

2 

5 = C/ono(r,t) = ^(l-^) (13) 

It is easy to see that the quantity w is a local velocity of the condensate. In fact, the total condensate wavefunction 
equals 



*o {r,t) = 7X0 exp 



(14) 



Taking into account (|^) and determining phase gradient, we arrive at expression (|T^) for the velocity. A special 
attention should be paid for the appearance of the Doppler vpx' term in Eq.(p^), p = {~ihV) being the operator of 
momentum. 

Note that, according to (|^), the density uq (r, t) is proportional to b^^ {t). Correspondingly, the correlation length 
of ^ ^ ?1q changes in time as ^ = (t) where <^o = ^ {2mUonQ (0, 0)) . The size of the system R (t) (see (|l0|)) 
changes in the same manner. Thus a ratio of these quantities ^ (t) /R{t) = S,o/Ro conserves for an arbitrary scale of 
evolution. 

The inequality S,o/Ro <C 1 is crucial for the Thomas- Fermi approximation. So, if this approximation is valid at 
the initial time, it remains valid for an arbitrary time. Since ^ (t) / R (t) <^ 1 , we have a quasi uniform problem at 
any time, considering excitations of the £, (t) < X <ti R {t) wavelengths. The scale of the uniform regions, according 
to Eq.([Tl|), changes with the typical time of t^ff ~ b{t) /b{t) . The well-defined excitations satisfy the inequality of 
ujteff ^ 1, i.e., satisfy the quasistationary condition. 

Let us consider the case of the fast transition from frequency to uji — ujq/ (3,(3 > 1. In this case the solution of 
Eq.(|) yields 

b^ it) = i (/52 + 1) _ i (^32 _ 2ujit. (15) 



Hence 



■^r^{f3^-l)sm2iu,t. (16) 



As we will see below, the excitations of A ^ give the main contribution into the damping. The energy of these 
excitations is of the order of UqUo (0, t) = fi/b^ (t) and, correspondingly, 

hujo (/32 - l)sin2wir 

It is clear that the ojteff ^ 1 condition holds if 2/i/?ia;o ^ /3. The last inequality can easily be fulfilled. Thus, we can 
seek for the solution of Eq. ( pT| ) within the quasiuniform and quasistationary approximation. 
Let us introduce a typical representation for operator x' 

x'/b = cxp {tkr^ aj:, (17) 

where aj: is the annihilation operator of a particle. We employ the known Bogoliubov transformation, e.g., |Q, for the 
operator 
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where bi and fej are the creation and annihilation operators of collective Bose excitations. As a result, the solution 



of Eq. ( 11 ) gives the excitation spectrum 
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ej:^hkv + ej:, e^^^^El-g\ i?- = 1^ + (19) 



The expression for the coefficients of the transformation is given by 



Let an external perturbation of a size of about d << Rq and localized near r = Tp be described by an effective 
interaction B{r — Tp). The Hamiltonian of the interaction of the gas with the perturbation has the form 

H,nt = j dPri>+{f,t)B{f-rp)i>{f,t) J cPrx+ {r,t)B{r - fp)x{r,t). 

Selecting terms corresponding to the single-particle excitations, we have 

1 /•.2, 



H^nt = J d'r[xQ{f,t)B{r~fp)x'{r,t) + hc]. (21) 
Let us transform this Hamiltonian, using (|l7|), (p^, and value of xo = b^/no 

H,nt = ^no{fp,t) % exp (zfcrp) ~ {k + (22) 

k 

If we take into account Eq.(|l9|), the probability per unit time of creating excitations equals {T — 0) 

W = ^no (fp, t) J \Bkf [uk - vkfh (sk + hkv^ . (23) 

In accordance with Eq. ( pO[ ) 
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2mefc 

Integration over the angles in Eq.(|2^) yields 

WiX-iiklfikv)) 



dipS (hkv cos if + ik) 



hkvij 1 — {ik/Tikvf' 

where {x) is the step function: {x) = 1 for a; > and {x) ~ for a; < 0. Then 

27rwm J £k L /-/*■; \2 

yl - {ek/fikv) 

As follows from this expression, creation of an excitation with the wavevector k occurs provided v > ik/hk, i.e., 
provided that the familiar Landau criterium holds. The generation of excitations is completely absent if v{fp,t) < 
c{fp,t). 

The local velocity of sound equals (see (p^)) 

g{r,t) _ Co / H ^ ' 



where cq = y^n/m. Taking into account the local magnitude of the gas velocity (|l^ ) and expression (^5|), we have 
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For 1 < ?7 < /?, the dependence of a on &^ {t) is a monotonically decreasing function with the inevitably existing region 
where the generation of excitations (a > 1) takes place. If 1 < 77 ^ /?, 

1 



If 77 < 1, the value a as a function of has a maximum at the point 

Calculating at this point, we can verify that the generation of excitations does not appear in the whole region of 
Vp < Rq \i (3 < v^- At the same time, the generation is absent at 7^< i?o/\/2 for an arbitrary value of /?. 

Suppose that Vp > Rq and let us calculate an integral in Eg. (p4|) . The values of k are limited by the condition 
which follows from the equality ik = Tikv, 

2m 



k<k,, k, = — \/i;2 _c2. (26) 



In the 2D case the effective vertex in the interaction Hamiltonian (|2^) at fc^d < 1 can be represented as 

5(0) 

Bk - y --^-r , (27) 

where b'^^ is the Fourier transform of the bare potential, d being a length of the order of the size of perturbation 
potential in the 2D case. If this value is large enough, we have 



~ mln(l/fcd) (^^^ 



i.e., it has only a weak logarithmic dependence on k. In the opposite case k^d > 1 Eqs. ( psD and ( Pq ) are not valid and 
the calculation of the transition amplitude requires a knowledge of the concrete form of the interaction potential. The 
most contribution into the integral ( p^ ) is given by the interval adjoined the upper limit k^. We can approximately 
take a factor of \Bk\^ out of the integral, putting fc* into the argument of the logarithm. Next, integral ( p^ can be 
calculated exactly as an integral over x — {ek/fivk) within the limits (c/v, 1) 

mno{fp,t)\B.^\^ 2 . c / c\ 
W = ^ 1 arcsm - 1 , B* = Bk ■ (29) 

Let us determine the total number of excitations created for a period of the condensate oscillations in the case of 
1 < 77 < /3. 

/ = 2 y dtW. (30) 
tl 

The time ti determines the moment when the boundary of the gas front reaches the point Tp. This value is found 
from the condition rp = b (ti) Rq. It follows from Eq.(f5|) that b^ {ti) « 1 + (woii)^ and ti m \/ "'f — ^/yjQ ^ The 



time ^2, found from the condition a = 1 (^sj), equals t2 ~ [7r/2 — 1/ (\/2?7)]/wi. Substituting Eq.(^ for the density 
into (^9|), one can see that the integral over time is dominated by the lower limit. From the physical point of view this 
is associated with the fact that the condensate density at point rp has a maximum during evolution just for t ^ ti. 
It follows from ( |2^ ) that c/v < 1/ (^y/2r]). Under this condition the factor in the square brackets in Eq.([2S|) is of the 
order of unity. After integration we find 

I ^ ^mno{0,0)\B,f —. (31) 



Substituting (28) into (^ij) and introducing the total number of particles N — i7ri?Q7io(0, 0) instead of no(0, 0), we 
find 
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W ln2(l/fc*d) ' i?o 

In the present case we have put huJo/fJ- <^ 1. This means that the number of excitations / and also the number of 
particles escaped from the condensate for the oscillation period are small compared with N. 

The total energy of the excitations generated during one oscillation period can be estimated as firj^I- Since the 
energy of the oscillating condensate is about fj,N, the relative energy loss for a period at 77 > 1 and /3 3> 1 approximately 
equals rf'I /N . Varying the initial condensate density (/i) or the 'impurity' position, one can change the damping rate. 

Treating a realistic problem in the quasi 2D case corresponding to the cylindric symmetry of the field configuration 
with the transverse parabolic potential and perturbation independent of the longitudinal coordinates, we arrive at 
the result of (^ ) with N equal to the total 3D number of particles. One can readily see that, if the condensate 
3D wavefunction is written as {x,y,z) — ^E'q {x,y)/y/L^ and normalized by the total number of particles, the 



-0 V-^i yT^J ~ ^0 

(21 

wavefunction is governed by the same Gross-Pitaevsky equation as in the purely 2D case with the replacement 
of C/q — > Uq"^^ — > U^^ /Lz- Since rig^' = Lz'n'^\ we have rS^'^U^^ = u'^^Uq^K The cylindric symmetry of perturbation 
results in creating collective excitations oi kz = alone. In addition, the problem of scattering proves to be purely 
two-dimensional and we arrive at Eqs.(|23) and (|32|). The quantities B^"^ and remain as 2D Fourier components. 



Note, that, if the laser beam parallel to the cylindric axis is used as a static perturbation (possibly, with the light 
sheets at the edges), an influence of the regions where the beam enters and leaves the condensate is unessential. 

In conclusion, we have investigated the generation of excitations in the oscillating condensate in a time-dependent 
parabolic trap at the presence of a static local perturbation. We have revealed the existence of two space regions for 
the position of an "impurity". In the first region of < Ro/\/2 the generation of excitations is absent at any scale of 
oscillations. In the second one of rp > Rq the generation is realized inevitably for an arbitrary set of parameters and 
we have found the generation rate. The results for the first region are a direct consequence of SF of condensate. For 
any position in the second region, there is a finite time interval when the local velocity of the condensate proves to 
be larger than the local velocity of sound, entailing violation of SF. Measurements of the damping of the condensate 
oscillations at different positions of the perturbation opens a possibility to study the manifestation of SF in the 
evolving condensate. 
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